arXiv: 1505.04513v2 [cond-mat.mes-hall] 21 Sep 2015 


Localization, delocalization, and topological phase transitions in the one-dimensional 

split-step quantum walk 

Tibor Rakovszky 

Institute of Physics, Eotvos University, Pazmany Peter setany 1/A, H-1117 Budapest, Hungary 

Janos K. Asboth 

Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, 

Hungarian Academy of Sciences, H-1525 Budapest P.O. Box 49, Hungary 

(Dated: Spring 2015) 

Quantum walks are promising for information processing tasks because on regular graphs they 
spread quadratically faster than random walks. Static disorder, however, can turn the tables: unlike 
random walks, quantum walks can suffer Anderson localization, with their wavefunction staying 
within a finite region even in the infinite time limit, with a probability exponentially close to one. It 
is thus important to understand when a quantum walk will be Anderson localized and when we can 
expect it to spread to infinity even in the presence of disorder. In this work we analyze the response 
of a one-dimensional quantum walk - the split-step walk - to different forms of static disorder. 

We find that introducing static, symmetry-preserving disorder in the parameters of the walk leads 
to Anderson localization. In the completely disordered limit, however, a delocalization transition 
occurs, and the walk spreads subdiffusively to infinity. Using an efficient numerical algorithm, we 
calculate the bulk topological invariants of the disordered walk, and find that the disorder-induced 
Anderson localization and delocalization transitions are governed by the topological phases of the 
quantum walk. 

PACS numbers: 03.67.Ac,73.20.Fz,03.65.Vf,05.60.Gg 


I. INTRODUCTION 

Discrete-time Quantum Walki 1 (or, simply, quantum 
walks) are quantum mechanical generalizations of the 
random walk. Their hallmark property is that on a reg¬ 
ular lattice they spread quadratically faster than ran¬ 
dom walks, i.e., ballistically rather than diffusively. This 
makes them valuable in quantum search algorithms! 2 , 
or even for general quantum computing-! Experiments 
on quantum walks range from realizations on trapped 
ionf®, cold atoms in optical lattice!^®, to light on an 
optical table®®, but there are many other experimental 
proposal J®®. 

The dynamics of a quantum walk is given by iterations 
of a unitary timestep operator, which can always be writ¬ 
ten in the form U = e~ lHeil , with H e s a hermitian op¬ 
erator. In this sense, a quantum walk is a stroboscopic 
simulator of an effective Hamiltonian This is a pow¬ 
erful theoretical concept that allows much of the physical 
intuition about lattice systems to be applied to quantum 
walks. As an example, consider quantum walks on reg¬ 
ular lattices: the maximum of the group velocity of the 
effective Hamiltonian translates directly to the velocity 
of ballistic expansion of the walk. 

In the presence of static (time-independent) disorder, 
quantum walks can lose their advantage over random 
walks in terms of the speed of spreading: they can un¬ 
dergo Anderson localization, whereby the mean-squared 
distance of the walker from the origin stays bounded even 
in the infinite-time limit. Besides theoretical 18 19 and 
numerical studies 2 ®, this effect has also been observed 
experimentally 21 . 


There are special cases where quantum walks can evade 
Anderson localization and spread indefinitely even in the 
presence of static disorder. Already the simplest one- 
dimensional quantum walk with angle disorder presents 
such a case: rather than being completely localized, it 
spreads subdiffusively 22 . This feature was explained in 
Ref.[22]by mapping the effective Hamiltonian of the quan¬ 
tum walk to chiral symmetric quantum wires (also see 
Ref. (221) ■ With an eye towards potential applications 
of quantum walks, it is important to understand under 
what conditions we should expect Anderson localization, 
and when delocalized behaviour, of disordered quantum 
walks. 

One of the key concepts that can help us understand 
when to expect Anderson localization in a quantum 
walkP, is that of topological phases. As noted in Ref. 1251 
the effective Hamiltonian H e g- of a quantum walk on a 
regular lattice can be engineered to be that of a topolog¬ 
ical band insulator®. If that happens, bulk-boundary 
correspondence® predicts that the quantum walk will 
host topologically protected edge states, whose number 
is given by a topological invariant of the bulk. These 
states can have a drastic influence on the time evolution 
of the walker if they have a large overlap with the initial 
state. 

Quantum walks have been shown to have a broader 
range of topological phases than that of their effective 
Hamiltonian- 2 ^ their topological invariants depend on 
details of how the timestep is performed. These invari¬ 
ants can be expressed as winding numbers of the bulk 
timestep operator over a part of the timesteji®®, or 
using a generalization of the scattering theory of topo- 
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logical phase®^. In this respect quantum walks are 
representative of the extreme limit of periodically driven 
systems-2UH1, a s is the (closely related) quantum kicked 
rotator 37 . In the case of the one-dimensional split-step 
quantum walk, topologically protected edge states not 
predicted by the invariants of the effective Hamiltonian 
have even been observed experimental!)®! the corre¬ 
sponding topological invariants have only recently been 
identified®. 

In this paper we explore the relation between Anderson 
localization and topological phases for one-dimensional 
split-step quantum walks. This is a broad family of 
quantum walks with chiral and particle-hole symmetries, 
which contains as a special case the simple quantum walk 
of Ref. [22. We use a cloning procedure to derive the 
real-space scattering matrix of the walk, which allows 
us to give simple and efficient formulas for the topologi¬ 
cal invariants as well as for the localization lengths. We 
find that uniform disorder in the rotation angles, which 
does not break the symmetries, leads to Anderson local¬ 
ization in the generic case. At maximal angle disorder, 
however, there is a disorder-induced delocalization tran¬ 
sition and the walk spreads (sub-)diffusively. We also ob¬ 
tain a simple interpretation of the delocalized behaviour 
of the simple quantum walk, by mapping it to a split- 
step walk at a boundary between topological phases. Fi¬ 
nally, we explore the effects of symmetry breaking dis¬ 
order using phase disorder, i.e., position-dependent but 
time-independent phase factor applied to the wave func¬ 
tion after every timestep. We find that when this disorder 
breaks the symmetries of the system it invariably leads 
to Anderson localization. 

This paper is structured as follows. In Sec. [TT] we re¬ 
mind the reader of the definition of the split-step quan¬ 
tum walk 2 ® 1 , of timeframes, and of the symmetries of this 
quantum walk. In Sec. [ITT] we derive our main result: the 
Lyapunov exponents of the split-step quantum walks, ob¬ 
tained through the real-space scattering matrix using a 
cloning procedure. In Sec. |IV| we apply this tool to treat 
uniform disorder in the rotation angles, find Anderson lo¬ 
calization or delocalization depending on the parameters, 
and give a new perspective on the delocalization of the 
simple quantum wall®]. In Sec.[v]we discuss two types of 
phase disorder: one which breaks both symmetries of the 
quantum walk and thus leads to Anderson localization, 
and one which breaks only particle-hole symmetry, and 
leads to the same type of behaviour as disorder in the 
rotation angles. In Sec. |VI| we draw some conclusions. 
We also include pedagogically important examples and 
calculations in Appendices. In App. [A] we discuss the 
pedagogical case of split-step quantum walks with binary 
disorder. In App. |B| we calculate the critical exponent of 
the localization-delocalization transition. In App. [C] we 
calculate the topological invariants of the split-step walk 
using the noncommutative generalization of the winding 
number'®. 


II. THE SPLIT-STEP QUANTUM WALK 

In this work we consider the split-step quantum walk 
on a one-dimensional chain. The state of the walker at 
each integer time t £ N is represented by a wavefunction, 

OO 

1^0))= E ip(x,s,t)\x,s). (1) 

X —— OO S=t4. 

Here x £ Z is the coordinate and s £ {f, 1} is the inter¬ 
nal degree of freedom, which we refer to as spin (in the 
quantum walk literature, the term “coin” is often used). 

The quantum walk consists of a sequence of three types 
of operations. Spin-dependent shift operations displace 
the walker, but do not mix the two spin components, 

Sf = X] (l a: + + |z4)(z4l); (2) 

X 

Si = E (|z,t)(z,tl + \x- 1, j)(z4|). (3) 

X 

Spin rotation operators rotate the spin about the y axis 
through a position-dependent angle, 

m = E E e- ie ^\x,s)(x,s\. (4) 

X S=t,4 

Finally, phase operators multiply the wavefunction by a 
position- and spin-dependent phase factor, 

P (0)=E E e~ i ^\x,s)(x,s\. (5) 

X s=t,i 

The quantum walk is defined by a short sequence of 
operations which is then periodically repeated. The effect 
of this sequence is represented by the unitary timestep 
operator, which for the split-step walk reads 

U(6 1 ,9 2 ) = S^R(0 2 )S t R(0 1 ). (6) 

The time evolution is given by 

m + l)) = U(0 1 ,0 2 )Mt)), (7) 

as represented in Fig. [T] 

The simple quantum walk, as, e.g., in Ref. 12~2[ has a 
single rotation operation per period. Its timestep reads 

u e ( 0 ) = S^Rtf). ( 8 ) 

This can be seen as a special case of the split-step quan¬ 
tum walk, with 02 = 0. 

We will also consider the phase disordered quantum 
walk. The timestep operator of this walk reads 

U p (<t>,6 1,0 2 ) = P(cl ) )S i R(8 2 )S t R(e i). (9) 
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FIG. 1: (Color online) A full timestep of the split-step walk, 
consisting of 4 operations: (1) rotation of the spin about the 
y axis through angle 20i, (2) displacement of the s = +1 
component of the wavefunction to the right, (3) second spin 
rotation through angle 202, (4) displacement of the s = — 1 
component of the wavefunction to the left. 


A. Timeframes, timestep operators, effective 
Hamiltonians 

There is a considerable freedom in specifying a quan¬ 
tum walk, i.e., a periodic sequence of operations, cor¬ 
responding to shifting the starting time of the period, 
which we refer to as changing the timefram^S. As an 
example, consider the operators 

u{e ll e 2 )' = R(e 1 /2)s i R(e 2 )s t R(e 1 /2) ] (io) 

U(6 1 ,9 2 )" = R(6 2 /2)S t R(9 1 )S i R(9 2 /2). (11) 

These both correspond to the split-step quantum walk 
as defined by Eq. only in different timeframes. 

Timestep operators describing the same quantum walk 
in different timeframes are related to each other by a 
unitary transformation. 

The effective Hamiltonian H c g of a quantum walk is 
defined as the logarithm of the unitary timestep operator, 

H eS = i\nU. (12) 

The branch cut in the logarithm is taken along the nega¬ 
tive real axis, and so all eigenvalues of the quasiener¬ 
gies E , are between — n and 7r. Since the unitary timestep 
operator U depends on the choice of timeframe (initial 
time of the period), the same quantum walk has many, 
unitary equivalent effective Hamiltonians associated with 
it. 


B. Symmetries and topological phases 

The split-step walk has both particle-hole symmetry 
represented by complex conjugation A', and chiral sym¬ 
metry, which places the system in Cartan class BDI 2 S. 


To see particle-hole symmetry of the quantum walk, 
notice that all matrix elements of the timestep operator 
U{9 1,02) are real (in position and <r z basis). Thus, 

KU(9 1 ,9 2 )K = U(9 1 ,9 2 ), (13) 


and therefore KH e ffI\ = — which is the defining 
relation of particle-hole symmetry. This holds for the 
split-step quantum walk in the two timeframes defined 
by Eqs. (10) and © as well. We remark that for a 
periodically driven particle-hole symmetric Hamiltonian, 
the symmetry is inherited by the effective Hamiltonian 
in all timeframeJ^]. 

To see chiral symmetry of a quantum walk explicitly, 
it is necessary to go to a chiral symmetric timeframe 28 . 
In the case of the split-step walk, there are two such 
timeframes, specified by Eqs. (10) and <©■ In these 
timeframes, we have 


<T X u(e 1 ,o 2 ya x = u(6 1 ,6 2 y\ 


(14) 


and, consequently, cr x H' eS a x = —which is the defin¬ 
ing relation of chiral symmetry. However, unlike particle- 
hole symmetry, the chiral symmetry requirement is non¬ 
local in timers. The addition of a phase shift operation 
to the timestep, as in Eq. ([9]) breaks both time-reversal 
and chiral symmetries. 

The presence of chiral symmetry enables us to assign 
bulk topological invariants to the quantum walk, simi¬ 
lar to those in static systems. Due to the periodicity of 
the quasienergy, however, edge states - eigenstates of the 
chiral symmetry operator - can exist at either E = 0 or 
E = 7r quasienergies. This means that there are two dif¬ 
ferent topological invariants, vq and u n , associated with 
the bulk. 

In the translation invariant case, when the angles 9\ 
and 9 2 are the same for all coordinates, we can obtain 
the the topological invariants by calculating the wind¬ 
ing number in both chiral timeframes given by (10) and 


(11). These two winding numbers can be combined to 
give the topological invariants and z^ 2 ^, resulting in 
the topological phase map shown in Fig. [2] 


III. LYAPUNOV EXPONENTS AND 
TOPOLOGICAL INVARIANTS BY CLONING 

In this paper we will be concerned with the effect of 
disorder on the topological invariants and edge states of 
the split-step walk. Since disorder breaks translation in¬ 
variance, the bulk topological invariants can no longer 
be obtained as fc-space winding numbers. There are two 
alternative approaches to the topological invariants for 
the di sorde red case: the one based on the scattering 
matrbPiEl and the one based on a reformulation of the 
winding number in real space, recently used for the dis¬ 
ordered SSH modeP^. In the following we detail the first 
approach. In App. 0 we briefly describe the second ap¬ 
proach, and compare the results obtained via these ap¬ 
proaches. 
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FIG. 2: (Color online) Phase map of the split-step walk, with 
gapped phases indexed by their pair of bulk topological in¬ 
variants, (vo,v n ). Along the continuous (dashed) lines, the 
quasienergy gap at E = 0 (E = n) closes. 


To define a scattering matrix, we have to use open 
boundary conditions, with two translationally invariant 
leads attached to the scattering regiorPS) To obtain leads 
that host the right number of propagating modes, we 
omit the rotations in semi-infinite parts of the system, so 
the timestep operator there simply reads U = S^S^, as 
shown in Fig. [3] 



U = S 4 S t U = Si.RaS^fRt U = S& 


Left lead Scattering region 


x = L + l 

-o—# 

Right lead 


FIG. 3: (Color online) To implement the scattering matrix for 
the quantum walk, we break up the periodic boundary condi¬ 
tions and attach two leads to the two ends of the system. The 
eigenstates in the leads are plane waves with linear dispersion 
relations. 

The scattering theory of topological invariants gives a 
simple way to write down the topological invariants of 
a bulk: they are related to the reflection amplitudes of 
the bulk at the quasienergies of the edge states. Chiral 
symmetry ensures that the reflection amplitudes at these 
energies are real and the topological invariants for a bulk 
are given bjP2 


v E = ^r(E). (E = 0 ,tt) (15) 

The scattering matrix of a quantum chain of L “slices” 
is usually obtained from the transfer matrix, which is 
the product of the L transfer matrices representing the 
effect of each slice. The product structure of the transfer 
matrix makes it straightforward to obtain the statistical 
properties of a system with uncorrelated disorder, if the 


transfer matrix of each slice is obtained from variables 
that are local to the site. 

There is a problem with applying the transfer matrix 
method directly to a split-step walk. During one timestep 
the walker makes excursions to nearest-neighbor sites and 
thus is not only affected by variables that are local to one 
site. 

Our approach to dealing with this problem involves an 
extension of the Hilbert space, by introducing N copies, 
or “clones” of the original split-step walk consisting of N 
substeps. We explain the details below for the N = 2 
case, the generalization is straightforward. 


A. Cloning 


The central idea of cloning, as shown in Fig. [4j is that 
we double the number of internal degrees of freedom at 
each site, thus creating two clones, and we break up the 
time evolution between these clones. The first half of 
the step takes place on the first clone and the second 
part on the second clone. At the beginning of each step, 
the walker is shifted from one clone to the other by an 
operator D. In formulas, 


C4loned = SRD 

(16a) 

© 

II 

(16b) 

R = /? , 0 R 2 

(16c) 


D = H E \x, a, l)(x , tr, 2| + \x , cr, 2)(:r, a, 1| (16d) 

X (T=t>I 


where the indices 1 and 2 refer to the first and second 
clone, respectively. This definition ensures that on each 
clone, (Ddoned) 2 is the same as the original timestep in 
some timeframe depending on which clone the walker 
starts from. 

In order to have the original time evolution, we need 
to start the particle from the second clone. In that case, 
each step in the cloned walk is equivalent to half a step 
in the original split-step walk defined in Eq. ([6]), 


U, 


cloned 


= u f 


splitstep* 


(17) 


The main advantage of cloning is that the walker can¬ 
not return to the site from which it started in a sin¬ 
gle step. This enables us to write down simple for¬ 
mulas for the stationary states. Note that according 
to Eq. ( [I7| ) the eigenstates of the original split step 
walk with quasienergy E will appear as states with 
quasienergy E/2 in the cloned walk. 

In general, a state in the cloned system has the form 


4>(x, s, l)|:r, s, 1)+i/j(x , s, 2)\x, s, 2) (18) 

x s =t4 

The VKIb 2) and ip(l, 1) components of the wavefunctions 
are unaffected by the shift operation S, while the ^(f 
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FIG. 4: (Color online) Cloning of the split-step walk. Start¬ 
ing from the second clone, the state of the walker after 2 1 
timesteps of the cloned walk is equivalent to the state of the 
original walker after t timesteps. 


for V£, 

vFS = O(£)0S T ^ + 4^t), (22) 

where we introduced the shorthand 

G x (e) = (e- fe l 2x2 - AS 0 )” 1 (23) 

for the resolvent of the matrix A® 0 . This enables us to 
relate O and O on adjacent sites. 


B. Real-space scattering and transfer matrix 

We arrive at what can be called the real-space scatter¬ 


ing matrix from Eq. (211, using Eq. (22), 


*l+i 


Tc t;E 


*t 


(24) 


with the components of the above matrix defined as 


t* = e i£ (4 t + 4 0 G*(e)Af); 

(25a) 

hr = e lt (4/ + A]?G x (e)A° ^ ; 

(25b) 

t' x =e is (A^ + Ai 0 G x (e)A° x l ); 

(25c) 

< = e is (A^ + At°G x (s)A^). 

(25d) 


, 1) component is shifted one site to the right and V(4- 
, 2) is shifted one site to the left (both these states are 
also shifted to the other clone). Thus we can divide the 


wavefunction into three parts, 


T—1 

A 

II 

(19a) 

4 = V’(a ; ,l,2), 

(19b) 

i>° = OO,i, 1 ), i>(x, t,2)). 

(19c) 


We then define the matrix A to be the part of the 
timestep operator preceding the shift in the basis where 
the wavefunction is ordered as Vx = (VJ , Vi , VS) : 

/ A n A n Af°\ / A t\ 

A = RD = A++ A-W- A±° = \A± , (20) 

yA°t A°+ A 00 / \A°J 


We can interpret Eq. (24) in the following way: the com¬ 
ponents and \Ed act as incoming ’’modes” that come 
to the site x during a specific time step from the left/right 
respectively, while an d 4+1 act as outgoing modes 
towards the left/right from the same site. This justifies 
the name real-space scattering matrix for the matrix ap¬ 
pearing in Eq. ( |24| . 

We can now obtain the reflection and transmission am¬ 
plitudes that relate out- to ingoing plane waves in the two 
leads. First, we combine the L individual 2x2 real-space 
scattering matrices of the sites by the usual combination 
rule of scattering matrices, and thus get the real-space 
scattering matrix of the whole scattering region. From 
the components ti,L and ti,£ of that matrix, the reflec¬ 
tion and transmission amplitudes are given by the fol¬ 
lowing formulas, 


where we introduced A mn that maps from sector n of 
the wavefunction to sector m, with m,n £ {f, i, 0}, and 
A m = (A m V A m V A m0 ). Each component depends on x. 

The equation for the components of the scattering 
states can be written as 


4+i4+i 

aw 



( 21 ) 


where e = E/2 is the quasienergy measured in the cloned 
system. The last component of Eq. (21) can be solved 


r(e) = e Je ti iL ; (26a) 

t{e) = e ie ti, L . (26b) 

The cloning of the system and the real-space scatter¬ 
ing matrices give us a useful method for calculating scat¬ 
tering amplitudes and thus topological invariants both 
numerically and analytically. This method can be easily 
generalized to more complicated quantum walk protocols 
involving multiple shift and rotation operators, or larger 
number of internal states. Later, in Sect. [V] we will use 
it for the phase disordered split-step walk as well. 
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Although for numerical work, the real-space scattering 
matrix is a practical tool because of its numerical stabil¬ 
ity, for analytical formulas, the real-space transfer matrix 
M x is more useful. It is defined by the relation 




£+1 

\]/i 


= 


^ X 

Vj/I . 


(27) 


For the split-step walk, the real-space transfer matrix 
at quasienergy E depends on the local angle parameters 
9j(x) through their sines and cosines, abbreviated as 

Sj = sin 0j(x); Cj = cos 6j(x), (28) 


with j = 1,2. The formula for the transfer matrix reads 


M X (E) = 


ClC 2 
a iE 


SlS 2 

—e iE / 2 si - e~ lE / 2 s 2 


—e _iB / 2 s 1 - e iE / 2 s 2 


—e 


-iE 


S\S 2 


(29) 


C. Lyapunov exponents, topological invariants, and 
localization length of the disordered split-step 
quantum walk 


To determine the topological invariants of the disor¬ 
dered split-step walk, we need the reflection amplitudes 
at quasienergies 0 and 7r, as per Eq. (151. At these 
quasienergies, the real-space transfer matrix of a single 
site, Eq. (29), is considerably simplified, 


M x ( 0) = e A * (o)CT * (30a) 

M x {n) = Ve x ^ ) ^V~ 1 , (30b) 


where V = V~ x = ( a x — <j y )/y/2 is a 2 x 2 unitary ma¬ 
trix. The parameters A^O) and A x (7r) are functions of 
the rotation angles # 1 , 2 ( 3 ;) of the site, 


= 1 (1 — sin# 1 (x))(l — sin# 2 (a;)) 

x{ ’ 2 (1 + sin# 1 (a;))(l + sin# 2 (a;)) ’ 

1 (sin#!(a;) - 1)(1 + sin# 2 (x)) 
x 2 (1+sin# 1 (a;))(sin# 2 (x) — 1) 


(31a) 

(31b) 


As seen from the Eqs. (30), at the quasienergies E = 


0,7r the real-space transfer matrices of the single sites 
all commute. Thus, the quantities A x (0) and Aa,(7r) are 
additive: A system of L consecutive sites is characterized 
by their sum, or equivalently, their average, 


A e = jJ 2 X ^ E ) (E = 0,n). 


(32) 


We will refer to this average A e as the Lyapunov expo¬ 
nent. 

The reflection and transmission amplitudes of the 
whole syst em a t the relevant quasienergies are obtained 
from Eqs. (26a) using the Lyapunov exponents as 


r{E) = — tanh(LA,E); (33a) 

t(E) = 1/ cosh(LA^) {E = 0,tt). (33b) 


In the thermodynamical limit of L —> 00 , if the Lya¬ 
punov exponents are nonzero, the effective Hamiltonian 
of the quantum walk is an insulator, and |f(E)| —> 1. In 
that case, the topological invariants are obtained from 
Eq. (pi), 


ve = -^sign(A s ). (34) 

For a translation invariant system, the Lyapunov expo¬ 
nent is the same as the additive parameter of a single site, 
A e = A X (E) as per Eq. ©. The topological invariants 
as a function of the global parameters #1 and # 2 , as read 
off from the plots of Fig. [5j agree with the previously 
known results of Fig. [2} 



FIG. 5: Color online) The negative inverse of the Lyapunov 
exponents, — 1/Ao and — 1/A^-, as functions of the rotation 
angles #1 and # 2 , as described in Eqs. ( |31[ ). Along the white 
lines, Xe 1, thus these points are on the topological phase 
boundary. Comparison with Fig.[2]shows that the signs of the 
Lyapunov exponents give the correct topological invariants as 
per Eq. p4|. 


While the signs of the Lyapunov exponents give us 
the two topological invariants, their absolute values tell 
us about the degree of localization of the states with 
quasienergies E = 0 and n. To see this we define the 
quasienergy dependent localization length: 


m 


2 L 

l n (\t(E)\ 2 /4) 


(35) 
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where it is assumed that the length of the scattering re¬ 
gion is L —► oo. From Eq. (33b) it follows, that for large 
L , the transmission amplitude at quasienergy E = 0 or 
7 r can be approximated by t(E) « 2e _i l AE L The local¬ 
ization length at these energies is therefore related to the 
Lyapunov exponent: 


m = 1/\X E \ (E = 0,n) 


(36) 


Transmission at quasienergies E = 0,7r decays expo¬ 
nentially with the system size whenever A® / 0. This 
shows that in any of the topological phases the quantum 
walk is insulating at both E = 0 and E = n. On the 
phase boundary however, where A e changes sign, £(E) 
diverges and the walk is no longer insulating at one of 
these energies. 

At E = 0 and n, the localization length £(1?) also 
defines the characteristic size of the edge states. To see 
this, suppose that we have an interface between two bulks 
characterized by different Lyapunov exponents: X r E for 
x > 1, and X l E for x < 0. We want to study the existence 
of zero energy edge states localized near this interface. 
If we chose the components of the wavefunction (4/j, \&q) 
at the interface to be an eigenstate of a x with eigenvalue 
±1, then according to Eq. (30), at L 1 distance from 
the interface, the wavefunction will be proportional to 
exp(±LA(j) in the right bulk and exp(=FLAo) in the left 
bulk, as shown in Fig. [6] Since £(0) = 1/Ao, this indeed 
means that the edge state decays into the bulk with the 
characteristic length £(0). Similar argument holds for 7r 
energy edge states with an appropriate choice of bound¬ 
ary conditions at the interface. The above argument also 
shows that the bulk-boundary correspondence holds for 
the topological invariants defined in Eq. (34), since the 


only way to create normalized edge states is if the sign 
of the Lyapunov exponents is different in the two bulks. 
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FIG. 6: Example for the wavefunction of a zero-quasienergy 
edge state near the interface of two disordered bulks with 
differing topological invariants (x < 0 : Ao > 0 and x > 
0 : Ao < 0). We chose (4'j,'F(;) to be (1,1), thus the two 
components of the spinor wavefunction (upper row, lower row) 
are real and positive, and the wavefunction is not normalized 
to 1. The components of the wavefunction simply acquire 
factors of e ± '^ B | as we move away from the interface, so that 
the edge state decays exponentially into both bulks. 


At a phase transition between different topological 
phases, at least one of the A e changes its sign. Thus, 
as we approach the phase boundary, the corresponding 


characteristic length scale diverges and the states with 
quasienergy E = () or n that were previously localized 
become delocalized throughout the whole system, i.e., 
their support will scale with the system size L , which in 
turn leads to high values of transmission probability at 
these energies. This delocalized behaviour leads to a sub- 
diffusive propagation of a walker started from the origin, 
as way previously noted in the case of the simple walk) 22 
In the limit of infinite system length, we can write 
down exact formulas for the Lyapunov exponents even 
in the disordered case. As already stated, by disorder 
we mean a probability measure /i( 0 1 , 62 ) given on the 
parameter space, according to which the angles 9\ and 
02 are chosen at each site. As we increase the length of 
the system, the sum ]T ' A X (E) performs a random walk, 
with the coordinate x playing the role of the time variable 
and X x (E) playing the role of the distance covered in the 
x-th step. The two topological phases correspond to the 
cases when the random walk of the Lyapunov exponents 
drifts to the plus/minus infinity, where A e can be thought 
of as a time-averaged drift velocity. For a system tuned 
to the topological phase boundary, the left and right drift 
terms cancel out and the random walk remains centered 
around the origin. Since the Lyapunov exponents are 
self-averaging, in the L —» 00 limit, the ’’time average” 
A e becomes the ensemble average, 


X E = j d/z(0r, 9 2 )X x (E,9 1 ,9 2 ) (37) 


According to Eq. (34), the sign of the above integral gives 


the topological invariant, while its absolute value defines 
the localization length as seen in Eq. (36) This means 


that Eq. (37), together with the definition of the Lya¬ 


punov exponents, Eq. (31), enables us to calculate the 


exact topological invariants and localization lengths for 
any disorder given in the form of /r(0i,02). This is the 
main result of this paper. 


IV. SPLIT-STEP WALK WITH UNIFORM 
DISORDER IN THE ROTATION ANGLES 

As an illustration of the ideas developed in the previ¬ 
ous Section we calculate the topological phase map of the 
split-step walk with uniform disorder. We take the rota¬ 
tion angles 9\ and $2 randomly and independently from 
a box distribution of width W centered around mean 
values (9 1 ), (0 2 ). The corresponding probability density 
function, 


Ml, (h) = (e« 0 i> + w- 9 1 )- 0(01 - (01> + W)) 

■ (0((0 2 ) + W - 0a) ■ 0(02 - (02> + w)), (38) 

is illustrated in Fig. [ 7 ] (For another illustrative example, 
that of binary disorder, see App. [A|) 

In Fig. [8j we show the effects of disorder on the zero- 
quasienergy Lyapunov exponent, more precisely, on the 












FIG. 7: (Color online) The split-step walk with uniform dis¬ 
order in the rotation angles. The parameters are taken from 
the yellow square of size 2W x 2W, with uniform probability. 
In Fig. [8] we move the middle of the distribution along the 
dashed line and vary the value of W from 0 to n. 



numerically integrate Eq. (37) to obtain the Lyapunov 
exponent Ao, and from it, the mean reflection ampli¬ 
tude —tanhAo as well as the localization length at 0 
quasienergy £o- We find that the boundary between the 
two topological phases, where Ao = 0 (dashed line in 
Fig. [8), is at (9\) = —vr/4, independent of the disorder 
W. Along this line, the system is critical, and the local¬ 
ization length diverges at quasienergy E = 0. 


A. Disorder-induced delocalization 

The system undergoes a disorder-induced delocaliza¬ 
tion transition in the strong disorder limit of W = n, 
where the rotation angles are taken from the whole pa¬ 
rameter space with equal probabilites. At this limit, 
which corresponds to the central point in Fig. [8j we find 
analytically Ao = A^ = 0. Thus, starting from a clean 
system, by increasing the strength of the disorder, we 
eventually reach a topological phase transition point, as 
shown in Fig. [9] This is accompanied by a delocalization 
of states with energies E = 0 and E = tt as described 
in the previous section. Thus, while at small values of 
the disorder strength, the quantum walk is Anderson lo¬ 
calized at all energies (we confirmed this separately by 
numerically calculating the localization lengths at other 
energies), strong disorder induces delocalization at the 
specific energies protected by the symmetries of the walk. 
This is similar to the case of the SSH modeP®. We note 
that the qualitative difference with respect to Ref. 1551 
where in the infinite-disorder limit a transition to a triv¬ 
ial localized phase was found, in we checked that using 
different disorder for 9\ and 9 2 leads to a transition from 
one topological phase to another, rather than delocaliza¬ 
tion. 


B. No Anderson localization on the critical line 


FIG. 8: (Color online) Topological phase map of the split-step 
walk with (do) = 7t/4 and uniform disorder in both rotation 
angles. The angle of the polar plot is the average (9 1 ) while 
disorder strength increases in the radial direction so that the 
middle of the circle corresponds to the strong disorder limit 
(see text), (a) Average reflection ampitude for a system con¬ 
sisting of a single site, — tanhAo. The dashed line separates 
the two topological phases with v 0 = — 1 and no = +1, the 
(red) dot indicating the values of the parameters used in Fig. 
|10| (b) Localization length, computed for a system of infinite 
size. Along the phase boundary, the localization length £(0) 
diverges. 


average reflection amplitude for a system of size 1, and 
on the localization length. For simplicity, here we set the 
mean second rotation angle to (9 2 ) = 7r/4, and we vary 
the mean of the first rotation angle, (9 1 ) from —7r/2 to 
7 t/2 (along the dashed line in Fig. [7|, as well as the dis¬ 
order strength W. At each point in this phase map we 


There are special cases for the split-step walk where 
disorder does not induce Anderson localization: if at 
some quasienergy the average Lyapunov exponent, de¬ 
fined by Eq. (37), vanishes. Uniform disorder in the ro¬ 
tation angles realizes such a special case if the quantum 
walk is on average on a phase boundary, i.e., if 


(9i) = ±(9 2 ) + mr (39) 

for some n £ Z, as seen in Fig. [T| In these cases, the local¬ 
ization length has to diverge at some critical quasienergy. 
As we show below, both the critical quasienergy and the 
shape of the divergence can be explained by a mapping 
to the disordered simple quantum walk. 

Obuse and KawakampS have found that the simple 
walk, defined in Eq. |8]), with uniform disorder in the 
rotation angle 9 , does not undergo Anderson localization, 
no matter how strong the disorder is. A key part of their 
explanation is that apart from the chiral and particle- 
hole symmetries, the simple walk posesses a sublattice 
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FIG. 9: (Color online) Topological phase change as a re¬ 
sult of increasing disorder, in a system of 200 sites, aver¬ 
aged over 1000 disorder realizations. Two branches of the 
dashed (red) curve: the highest/lowest quasienergy in the 
lower/upper band. The distance between the branches - the 
gap around zero quasienergy - gradually closes as disorder 
increases. Solid (blue) curve: the average zero quasienergy 
reflection amplitude. As the disorder increases, this remains 
quantized at r( 0) = 1 as necessitated by Eq. (151 long after 
the gap closes. As the disorder strength is increased even 
further, (r(0)} begins to change and eventually at W = n 
it reaches 0, signifying a phase transition point. The (blue) 
points represent r( 0 ) in the 1000 individual realizations whose 
average is the solid blue curve. 


symmetry, 


AUA = ~U A= \ x )( x \~ k)(a;|, (40) 

x even x odd 

because the walker can only hop from even to odd sites 
or vice versa in one step. Using this extra symmetry, 
they argued (for more details, see Ref. [23j) , that Anderson 
localization is avoided because the localization length £ of 
the simple quantum walk diverges at E = ± 7 r/ 2 , scaling 
as 


£(U) = £ 0 | log(&Er)|, 


(41) 


where SE = E — E CI i t is the distance from the critical 
quasienergy Ucrit = ±7 t/ 2, and r is the mean free time. 
As illustrated in Fig. [lOj the same effect occurs in split- 
step walks, if Eq. ( |39| ) is fulfilled. We show this explicitly 
in the rest of this section. 

Using the sublattice symmetry property of Eq. (401, a 
split-step quantum walk can be understood as a doubled 
sequence of the simple quantum walk. Consider 



FIG. 10: (Color online) Divergence of the localizaton length 
£( 1 ?) as a function of quasienergy A at a critical point with 
mean angles (0}i = —( 6)2 = — 7 r /4 and disorder W = 7 t /2 
(see Fig. | 8 ]/(b)). (a) The average localization length diverges 
at E = 0. (b) The scaling of the localization length near the 
critical energy. The dotted line indicates the curve (41 1 fitted 
to the numerical data, which gives r = 0.216 and £0 = 0.322. 
We used a system of 100 and 400 sites for figures (a) and 
(b) respectively. In both cases we averaged over 100 disorder 
realizations. 


This mapping is closely related to that used in Ref. 133) to 
realize a split-step quantum walk as a periodically driven 
Hamiltonian. 

The mapping of Eq. (42) shows that if (0i) = ( 62 ), the 
correlation length of the split-step walk has to diverge at 
quasienergy E crit = n following the scaling of Eq. (Ell). 


By Eq. (42), two timesteps of a simple quantum walk 
with angle 0 can be understood as a single timestep of a 
split-step quantum walk with 9\ = 62 = 8 . The introduc¬ 
tion of uncorrelated box disorder in the rotation angle 9 
of the simple quantum walk translates to uniform disor¬ 
der in the angles 9\ and 62 of the quantum walk, with 
(9i) = ( 02 )- The doubling of the timestep moves the 
critical quasienergy from ±7r/2 to ± 7 r. 

A mapping between different split-step quantum walks 
shows that if ( 0 i) = ( 6 * 2 ) ± vr, the divergence of the 
localization length follows the same functional form of 
Eq. (41), with E cr jt = 0. Changing 9j{x) —> 9j(x) ± 7 r 
for all x, for either j = 1 or j = 2 , results in an overall 
factor of —1 for the timestep operator U. This shifts the 
divergence of the localization length of the quantum walk 
by 7 r, to E clit = 0. 

A mapping between the transfer matrices of different 
split-step quantum walks shows that if (9 1 ) = — ( 02 )+wr, 
the localization length diverges according to Eq. (41), 
with shifted E cr i t . According to Eq. 


changing either 
9\(x) to — 9i(x) or 0 2 (x) to —9 2 {x) changes the transfer 
matrix M x at x to M' , with 


(x , ,s , \[S i R(9 2 )S t R(9 1 )} t \x,s) = 

(2a/, s , |[S , i?(0 , )S’R(0 , )] t |2z, s), (42) 

with the position dependent rotation angles 0 '{x) defined 
by the relations 


(43) 

(44) 


M' X (E + n) = - (j J) M X {E) ( J) . (45) 

The transfer matrix M' X (E + n) is a unitary transform 
of M X (E), up to the unimportant factor of —1. Since 
the transformation is independent of x, the Lyapunov 
exponents at E of the original walk are equal to those of 
the transformed walk at E + ir. This shows that if (9i) = 
— ( 02 ), the localization length diverges at E CI i t = 0, while 


9'{2x) = 6i(x); 
9'(2x — 1) = 02(a;). 
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if (9 1 ) = — ( 62 ) ±7r, it diverges at E CI j t = 7r, following the 
scaling of Eq. (41). 

Our results also shed new light on the absence of An¬ 
derson localization in the simple quantum walli 22 ! As we 
have shown, the simple quantum walk with disorder in 
the rotation angles is not Anderson localized because it 
constitutes a disordered split-step quantum walk tuned 
to a topological phase transition point. 


V. SPLIT-STEP WALK WITH PHASE 
DISORDER 

We now consider the phase disordered split-step walk, 
introduced in Eq. <©• Since phase disorder breaks both 
chiral and particle-hole symmetry, we expect that it in¬ 
duces Anderson localization at all quasienergies. There 
is a way, however, to add phase disorder to the split-step 
walk and keep chiral symmetry: in that case, we expect 
to see localization-delocalization transitions as with angle 
disorder in the previous Section. We discuss both types 
of phase disorder, and illustrate our results by numerical 
examples that can be compared directly with those on 
angle disorder of the previous Section. 

The simplest way to introduce phase disorder is to mul¬ 
tiply the wavefunction of the walker at the end of each 
timestep by a position- and spin-dependent phase factor 
<j>(x, s), chosen randomly and independently at each site, 
as defined in Eq. ®>- For the examples in this section we 
used an extra restriction of 0(cc,t) = — <j)(x,-l), whereby 
the phase operator reads 

P(<t>) = J2 \x)(x\®e~ i,t ‘( x) ' T *, (46) 


with the phase chosen from an interval [—A<p, A</>] with 
uniform distribution. Just as with the more general phase 
disorder, due to this extra operation, both particle-hole 
symmetry and chiral symmetry of the quantum walk are 
broken. Thus, in the presence of phase disorder, there 
are no localization-delocalization transitions. 

To highlight the role of chiral symmetry, we also con¬ 
sider adding phase disorder to the split-step walk in a 
chiral symmetric way. This requires two phase operators 
per timestep, 

01,0 2 ) = P{ ( j ) )S i R(e 2 )S t P((t>)R(e 1 ). (47) 


The second phase operation restores chiral symmetry 
since a x P{(j))a x = P(0) -1 . Repeating the calculation 
of the real-space transfer matrix, Eq. (29), with R(9\) 
replaced by P((j))R(9i)P(<j>), we find that this matrix is 
just multiplied by a factor of e 2l ^ x \ For example at zero 
quasienergy, 


M x (0) = e 2 ^ (x W 0 > CT *, (48) 

and similarly for E = n. The extra phase factor drops out 
from both the reflection amplitude and the localization 


length so that the description we gave in the previous 
section is unaffected. 


To show the effects of phase disorder numerically, we 
have calculated the localization lengths £ for a range of 
angle disorders and phase disorders, shown in Fig. [lT 


For easy comparison with the results of the previous Sec¬ 
tion, Fig. [8j we fix (0 2 ) = 7t/4, and show the localization 
length £(0) as a function of the mean first rotation an¬ 
gle (9 1 ) and of disorder, introduced in equal measures 
to both rotation angles and to the phases (f>, setting 
W = A </>. Thus, the perimeter of the plot of Fig. El 
with W = 0, corresponds to the perimeter of the plots of 
Fig. IJr). As expected, symmetry breaking phase disor¬ 
der, Fig.[lTJi), destroys the delocalization transition and 
for large values of A <fi, the states with 0 quasienergy are 
localized for all parameters. Phase disorder that respects 
chiral symmetry, Fig. E b ) , however, has the numeri¬ 
cally obtained localization lengths taking large values at 
( 9 i ) = — 7t/4, which is the phase boundary between topo¬ 
logical phases, and where the theory (shown in Fig. [8]) 
predicts divergences. 



(a) (b) 


FIG. 11: (Color online) Average localization lengths at 
quasienergy E = 0 for phase disordered quantum walks 
with equal disorder strengths for the rotation angles and 
the random phases, Arf> = W. (a) Single phase shift. The 
timestep operator is U = P(<f>)SiR( 92 )SfR( 9 i) thus chiral 
symmetry is broken, (b) Double phase shift, with timestep 
U' = P(4>)SiR(9 2 )S-fP((j>)R(di), exhibiting chiral symmetry. 
Single phase disorder leads to topological triviality and causes 
all states to be localized while double phase disorder restores 
the phase map seen in Fig. [7] with a delocalization transition 
at the topological phase border. The localization lengths were 
calculated using a system of 80 sites and averaged over 100 
realizations. 


We also show the effects of symmetry breaking and of 
chiral symmetric phase disorder on the time evolution of 
the split-step quantum walk directly in Fig. [T2| Here we 
show the position variance of a particle started from the 
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origin after t timesteps, 

var(x) = (^(t)\x 1 2 \i/j(t)) - \{ip(t)\x\il){t))} 2 . (49) 

We used the strong disorder limit in the rotation angles, 
W = 7r, and used phase disorder with amplitude A <j> = 
7r/4. For a single phase shift (symmetry breaking phase 
disorder), the walker is localized and var(x) —> constant 
for large t. However, in the double phase shift case, when 
chiral symmetry is restored, we see the subdiffusive be- 
vaviour discussed in the previous Section, with the walker 
slowly spreading through the system. These results are 
in accordance with those from the localization lengths. 



FIG. 12: (Color online) Variance of the position of a quan¬ 
tum walker started from a single site as a function of time 
for single and double phase disorder, averaged over 10000 dis¬ 
order realizations, with angle disorder strength W = n and 
phase disorder strength A <f> = 7r/4. For a timestep containing 
a single phase shift (green thin curve), the walker is localized 
and the variance tends to a constant, while for the chiral sym¬ 
metric case defined in Eq. ( |47| ) (red thick curve), the walker 
spreads subdiffusively. In both cases, we fitted the function 
var(x) = a ■ t b to the data points with t > 1000 and found 
that the exponent b is 9, 57 • 10 -4 in the first and 0,107 in the 
second case. The dashed lines show the fitted curves. 


VI. CONCLUSIONS 

We considered the effects of disorder on the localiza¬ 
tion properties and on topological phases of the one¬ 
dimensional split-step quantum walk. We introduced 
an effective numerical tool (the cloning trick), which al¬ 
lowed us to calculate the scattering amplitudes for the 
split-step walk and efficiently calculate the topological 


invariants proposed in Ref. j32J We then showed theoret¬ 
ically and investigated numerically various localization- 
delocalization transitions that occur whenever this sys¬ 
tem is tuned to a critical point at a topological phase 
transition. We have shown using a mapping that the 
subdiffusive spreading of the simple quantum walk with 
angle disorder 22 can be understood in this framework. 
We have shown that angle disorder generically localizes 
the split-step quantum walk, but that complete disorder 
in the rotation angles places it in a critical state with 
subdiffusive instead of localized dynamics. Finally, we 
illustrated the importance of symmetries on localization- 
delocalization through the example of phase disorder. 

It is interesting to compare our results on the one- 
dimensional split-step quantum walk (ID) with those 
obtained in Ref. [24] regarding the two-dimensional split- 
step quantum walk (2D). In the 2D case, the topological 
phases did not require any symmetry of the system, and 
so phase disorder did not destroy the topological phase. 
Similarly to angle disorder in the ID case, phase disorder 
in the 2D case was found to lead to Anderson localiza¬ 
tion, except when the system was tuned to criticality 
(as in the case of the Hadamard walk). The disorder- 
induced delocalization transition however was reached in 
both the ID and the 2D case by using complete disor¬ 
der in the rotation angles. In the 2D case, it was found 
that angle disorder alone does not lead to Anderson local¬ 
ization, possibly related to the presence of particle-hole 
symmetry. The effect of particle-hole symmetric disorder 
remains to be studied in the ID case. 

Our results show how the understanding of topologi¬ 
cal phases of quantum walks can help interpret their be¬ 
haviour under different types of disorder. This could be 
important to identify which types of quantum walks are 
practical for information processing purposes, and which 
types of disorder it is crucial to supress in such applica¬ 
tions. 
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Appendix A: Split-step walk with binary disorder 

In the main text we used uniform disorder as an il¬ 
lustration of the general formalism. In this section we 
present another simple example that can be treated an¬ 
alytically and is useful to obtain intuition regarding the 
disordered quantum walk. Consider a split-step quantum 
walk where the rotation angles 9\ and 62 can take two 



13 


different values: Of, Of or Of , Of - At each site, we chose 
one of these two set of values, with probabilities q and 
1 — q, so that the corresponding probability measure is 


The critical line <? C rit(?7i) for E = tt is shown in Fig. 14 


along with the numerically calculated average reflection 
amplitudes at quasienergy n. 


v( 0 i, 0 2 ) = q 6 ( 0 1 - 0 f) 6 ( 0 2 - 0 f) 

+ (1 -q)S(0 1 -0f)6(0 2 -0f). (Al) 
We fix the point A in parameter space as 

Of = 0.625tt; Of = -0.125tt, (A2) 

and choose point B to lie on a straight line that is parallel 
to one of the phase borders and goes through point A , 

Of = Of — m n; Of = Of + m it, (A3) 

where the parameter to € [0,1] measures the distance 
between points A and B in parameter space as shown 
in Fig. [13] When m = 0, the two points coincide while 
to = 1 corresponds to the case when their distance is 
maximal. 



FIG. 13: (Color online) Split-step walk with binary disorder. 
The red and blue dots represent points A and B in the case 
to = 1 , while the dashed line shows the possible values of B , 
while m goes from 0 to 1. 

At to = 3/8, the point B crosses the line where the 
gap at E = tt closes and the invariant changes. Thus, 
for values to > 3/8 the two limits q = 0 and q = 1 belong 
to two different topological phases. We want to find the 
exact value g cr jt for each to where the phase transition 
occurs. Suppose that we have a number of La sites with 
parameters 0 f 2 in our system and a number of Lb sites 
from the point B. The condition for the phase transition 
can then by written as 

Xe = LaX^ + Lb Xf = 0. (A4) 

From this we acquire the critical value of the mixing prob¬ 
ability q , 

1 1 , AC , 
q ciit 1 1 Lb. 1 _ A a ' ( A5 ) 


(r(vr)) 



Distance in parameter space m 


FIG. 14: (Color online) Topological phases of the split-step 
walk with binary disorder. Color coding shows the reflection 
amplitudes at E = n for a system of 40 sites averaged over 
100 samples. The green dashed line s hows the critical mixing 
probability q CI it calculated using Eq. (A51. The peak at to = 
5/8 corresponds the the case when the point B is in the middle 
of a topological phases and diverges. 


This binary disordered model, while somewhat un¬ 
physical, shows the importance of the Lyapunov expo¬ 
nents defined in Eqs. ( |30| ). They serve as weight factors, 
determining how much a given site contributes to the 
overall topoogical invariants of the whole system. 


Appendix B: Critical exponent for uniform disorder 


In this section we will give the critical exponent with 
which the localization length £(0) or £(rr) diverges when 
we approach the topological phase boundary in parame¬ 
ter space. There is also a different critical exponent when 
we are at a phase boundary: then the function £(E) di¬ 
verges as the quasienergy approaches 0 (or tt). This was 
discussed in section m of the main text. 

Let us first look at the clean case. Near the phase 
boundary 0 \ = —02 we can expand the function Ao in 
the small parameter 6 = 0i + 0 2 . From the definition 
Eq. (31) we get that Ao oc 8 so that £(0) oc 5~ l . This 
means that in the clean system the critical exponent is 
v = - 1 . 

When uniform disorder is applied, we can approach 
the phase boundary in two different ways as shown in 
Fig. [8] either by changing the average values the rotation 
angles (similarly to the clean case) or by increasing the 
disorder strength until we reach the strong disorder limit 
as described in the main text. Numerical evaluation of 
the integral (37) in the vicinity of the phase boundary 
verifies that the critical exponent remains u = — 1 in 


both cases as seen in Fig. 15 
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FIG. 15: Behavior of the localization length £(0) near the 
topological phase boundary, (a) Split-step walk with uniform 
disorder, disorder strength set to W = 0.17T. The average 
rotation angles approach the gap closing line 9i + 62 = 0. (b) 
Near the strong disorder limit, W ~ n. In both cases, we 
found a critical exponent v = — 1. 


Appendix C: Comparison with the real-space 
winding number method 


As we mentioned in the main text, there is an alter¬ 
native way to calculate topological invariants for a disor¬ 
dered system, apart from the scattering matrix approach 
that we used to study the split-step walk. This other ap¬ 
proach, based on non-commutative geometry tools, was 
used to study two dimensional topological insulators at 
strong disorder, and, more recently, the disordered SSH 
modePSl. In this formalism, the winding number of a one 
dimensional topological insulator with only two bands 
can be calculated by the expression 

u = -Tr(Q_ + [X,Q + _])/L, (Cl) 


where X is the position operator, Q _= r_QT + and 

Q _| = r + QT_, with Q being the flat band version of 

the Hamiltonian H , defined by replacing each eigenvalue 
of H with its sign: Q = sign (H). The operators F + and 
T_ are the projectors associated to the chiral symmetry 
operator T by T = T + — T_. 

Evaluation of the winding number, Eq. (Cl), for a fi¬ 
nite size system with periodic boundary conditions re¬ 
quires non-trivial approximations as described in Pro- 
dan’s papei®0. In this case, the resolution in quasi- 
momentum space for a system of L sites is given by 
A = 2^t/L. We approximate the differential d k A{k ) for 
some function A{k) of the quasimomenta with discrete 
differences in the following way: 


<? 

dkA(k) -A 5 k A(k n ) = ^2 Cm[A(k n +mA)—A(k n —mA)], 

m—1 

(C2) 

where q is of order L/2. 

If we assume that A{k) can be expressed as a Fourier 
series, then it is enough to find the coefficients c m that 
give a good approximation for functions of the form e lkx . 
For these, the formula ( |C2[ ) gives 

(d k - 5 k )e ikx = i(x - 2^c TO sin(roAz))e ifcx . (C3) 

m 

We can make the above expression disappear in the first 
2q — 1 orders of Ax by choosing c m to be the solutions 


of the equation 


Mc T = —(1,0,0,..., Of; 


M nm = m 2n 1 

n, m = 1, 

• •• ,Q • (C4) 

In this case, ( d k — 5 k )e xkx = 

0{A 2q ~ l ). 

By choosing 

q = L/2 we can make the error 

of the approximation to 


be of order £>(A i_1 ). 

For the split-step walk, 
both chiral symmetric timeframes defined in Eqs. (10) 
and (11). This gives us two winding numbers v' and v" 


Eq. (Cl) can be used in 


The top olog ical invariants are given as combinations of 
these twcPS) 


^0 




u' + v" 
2 


2 


(C5a) 

(C5b) 


These equations, along with Eq. (Cl) enable us to cal¬ 


culate the topological invariants numerically. Below, in 


Fig. 16 we show the results for z/o, along with those per¬ 
formed using the scattering matrices described in the 
main text. The two methods yield the same results qual¬ 
itatively. However, we note that while the reflection am¬ 
plitude of the finite system remains close to the value it 
has in the thermodynamic limit even for relatively strong 
disorder, the real-space winding number is more easily 
affected by numerucal inaccuracies resulting from small 
system sites. The scattering matrix method is also more 
efficient numerically (it scales linearly with the system 
size). 


o 



FIG. 16: (Color online) For a disordered quantum walk, the 
bulk winding number (solid red) and the average reflection 
amplitude (dashed blue) give similar predictions about the 
topological phase. We used a system of 100 sites and aver¬ 
aged over 1000 disorder realizations. For the bulk winding 
number, Eqs. (Cll were used, while the reflection amplitude 
was obtained by iterating the scattering matrix. 
















